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Abstract 

We propose a new method for biack hoie spin measurement. In this method, we consider 
a gas biob or ring faiiing onto a biack hoie from the marginaiiy stabie orbit, keeping its initiai 
orbitai anguiar momentum. We caicuiate the gas motion and photon trajectories in the Kerr 
space-time and, assuming that the gas biob or ring emits monochromatic radiation, carefuiiy 
examine how it is observed by a distant observer. The iight curve of the orbiting gas biob is 
composed of many peaks because of periodic enhancement of the flux due to the gravitational 
lensing and beaming effects. Further, the intensity of each peak first gradually increases with 
time due to the focusing effect around the photon circular orbit and then rapidly decreases due 
to the gravitational redshift, as the gas blob approaches the event horizon. The light curve 
of the gas ring is equivalent to a superposition of those of the blobs with various initial orbital 
phases, and so it is continuous and with no peaks. The flux first gradually increases and then 
rapidly decays, as in the blob model. The flux variation timescale depends on the black hole 
spin and is independent from the inclination angle, while time averaged frequency shift have 
dependences of both effects. We can thus, in principle, determine spin and inclination angle 
from observations. The observational implications and future issues are briefly discussed. 
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1 Introduction 

Observational proof of the blaek hole spaee-time is one of the most outstanding issues in physies and 
astrophysies. This ean lead to a proof of the existenee of the event horizon and to a eritieal observa¬ 
tional test of the general relativity theory proposed 100 years ago. It is known in general relativity that 
blaek hole spaee-time is uniquely deseribed by a blaek hole mass, M, and a spin parameter, a (other 
than eharges, whieh are never important in the astrophysieal eontext). Masses ean be relatively easily 
estimated by observing the motions of stars or gas (Shahbaz et al. 1999; Ghez et al. 2005; Orosz et al. 
2011). The observed targets may not be neeessarily elose to the blaek hole, as long as they are bound 
to its gravitational potential. Blaek hole spins are, on the other hand, not easy to measure, sinee they 
only influenee the spaee-time near to the blaek hole so that full eonsiderations of general relativistie 
effeets should be neeessary. 

Until now, several methods have been proposed for the spin measurements. There are three 
major methods: (i) eontinuum speetrum method (Hanawa 1989; Li et al. 2005; MeClintoek et al. 
2011; MeClintoek et al. 2014), (ii) line speetrum method (Kojima 1991; Laor 1991; Tanaka et al. 
1995; Duro et al. 2011; Steiner et al. 2011; Reynolds 2014), and (iii) quasi-periodie oseillation 
(QPO) method (Kato 2001; Strohmayer 2001; Rezzolla et al. 2003; Remillard 2005; Kato et al. 
2008). 

We wish to note, however, that all of these methods have uneertainties and/or are based on 
the eritieal assumptions that are not always easy to prove. The eontinuum speetrum method (i), 
for example, eritieally depends on the inelination angle. It also has an uneertainty in the speetral 
hardening faetor in the emission from an aceretion disk. The line speetrum method (ii) is sensitive 
to the emissivity profile assumed for the illuminated flux, as well as the iron ionization fraetion as 
a funetion of radius. Further it is not obvious how to separate line speetra from eontinuum ones. 
Moreover, it is assumed in both of methods (i) and (ii) are eonstrueted on the assumption that the 
inner edge of the disk is at the radius of marginally stable orbit, but this remains as a eontroversial 
issue (see Shidatsu et al. 2014). In the QPO method (iii) no widely aeeepted model is available for 
the origins of quasi periodie oseillations. Further, several different modes of QPOs are known and it 
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is difficult to identify. 

We should also point that spins estimated by different methods are sometimes ineonsistent. 
For example, the spin of GRO 1655-40 estimated by methods (i) and (iii) are 0.65 < a/M < 0.75 and 
0.31 < a/M < 0.42, respeetively (Shafee et al. 2006; Kato et al. 2008; MeClintoek et al. 2011). 
The spin of GRS 1915-1-105 is estimated to be a/M > 0.98 [method (i)], 0.54 < a/M < 0.58 or 
0.97 < a/M < 0.99 [(ii)], and a/M < 0.44 [(iii)], respeetively (Kato et al. 2008; Blum et al. 2009; 
MeClintoek et al. 2011). In short, all the existing methods are far from being eomplete. Therefore, it 
is required to eonstruet new and independent methods for spin measurement to eheek other methods. 

That is, we eonsider the gas blob or ring falling into the blaek hole from the marginally stable 
orbit (or the inner edge of an aeeretion disk) with a finite angular momentum and earefully examine 
how the gas blob or ring is observed by a distant observer. Radiation from the gas is known to undergo 
five relativistie effeets whieh are summarized in table 1 (see Karas et al. 1992; Dovciak et al. 2004a). 
We ealeulate the gas motion and photon trajeetories in the Kerr spaee-time by general relativistie 
ray-traeing method to seek for good observational indieators of the blaek hole spin. Note that our 
approaeh is distinet from that by Dovciak et al. (2004b) who were eoneerned with radiation from an 
aeeretion disk whose inner radius ean be set below the marginally stable orbit, whereas we eonsider 
the gas blob or ring whieh is separated from an aeeretion disk and falls from the radius to the blaek 


hole. 


The plan of this paper is as follows: in seetion 2, we introduee the basie equations to deseribe 
gas motions and ray trajeetories, explain our models for the gas blob and ring, and deseribe methods 


of ealeulations. In seetion 3, we show numerieal results and explain the physies underlying the key 
features of the light variations. On the basis of these simulation results we propose in seetion 4 our 
method for determining a spin a and inelination angle i. Seetion 5 is devoted to diseussion on the 


distinetive features of our method, its observational implieations, and future issues. 

2 Models and methods of numerical calculations 

In this seetion, we first introduee the basie equations (subseetion 2.1). Next, we explain our models 
of an infalling gas blob (subseetion 2.2) and ring (subseetion 2.3). Finally we deseribe methods of 
numerieal ealeulations (subseetion 2.4). 

2.1 Geodesic equations of gas particies and photons 
In Boyer-Lindquist eoordinates, the blaek hole spaee-time is 
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Table 1. Five relativistic effects on radiation emitted by gas near a black hole. 



Relativistic effect 

Physical cause 

Remark (reference) 

(I) 

Gravitational lensing 

Focusing by light bending 

towards a gravitating source 

(Karas et al. 1992) 

(II) 

Beaming 

Energy boost by relativistic 

motion of particles 

(Karas et al. 1992) 

(III) 

photon circular orbit 

Focusing around the 

photon circular orbit, rph 

key feature in the early stage 

(IV) 

Focusing by the 

frame dragging 

Focusing by frame dragging 

due to black hole spin 

effective when a/M is high 

(V) 

Gravitational redshift 

Energy loss by light travel 

from deep potential well 

key feature in the late stage 


/ 2 2 2Ma^rsin^0\ . 2 / 17-2 

+ H--- sin^ Odcj), 


( 1 ) 


( 2 ) 


with 

A = — 2Mr + a^, 

E = + a^cos^O, 

where we shall use “geometrized units”, in whieh the gravitational eonstant, G, and the speed of 
light, c, are set to be unity. The spin parameter of the blaek hole is a = J/M, where J is the angular 
momentum of the blaek hole. 

Around the blaek hole, motions of a partiele (with mass of /i) and photon (with no mass) are 
determined by the geodesie equations (Bardeen et al. 1972): 

E— = —aiaEsm. 9 — L) + - ---, 

dX A 




+ 


aT 

A“’ 


( 3 ) 


with 


E = —pt = const, 
L = p^ = const, 


Q =Pe + cos^ 6 




P<p 

sin^^J 


= const. 















(4) 


T = E{r^ + a^) - La, 

Vr = T‘^- + (L - aEf + Q], 



where is the 4-momentum of a test partiele, and A is related to the proper time of the test partiele, 
r, by A = r/yU. Note that A is the affine parameter in the ease /i —?■ 0. In addition, E, L, and Q are the 
energy, the angular momentum, and the Carter eonstant of a test partiele, respeetively (Carter 1968). 

As for a photon trajeetory, we set /i = 0 in equations (3) and (4). 

2.2 Gas blob model 

We postulate the situation that a part of the innermost region of an aeeretion disk is stripped off, 
thus forming a gas blob, and that it starts to fall onto the blaek hole, keeping its original angular 
momentum. We assume that the gas blob has the following properties: 

1. It has a spherieal shape and its eharacteristic radius is i?biob> which has no time dependence. 

2. We neglect self-gravity of the blob. 

3. The gas blob emits radiation with monochromatic frequency, pq, in its inertial frame. This assump¬ 
tion is valid even in the case of multicolorwavelength radiation, such as the blackbody radiation, if 
the frequency width of the spectrum is sufficiently small. Its emissivity decreases with an increase 
of a distance from a center of the blob, obeying the Gaussian function. 

4. We neglect self absorption within the blob. 

5. The initial position of the center of the gas blob is at (?",(/)) = (0.98rms,0). Here r^s is the radius 
of the marginally stable orbit. 

6. The motion of the center of the blob follows the one particle orbit given by equations (3) and 
(4). The blob falls to the black hole on the equatorial plane, keeping the constant energy, E^s, 
and angular momentum, L^s, where i^ms and Lms are the energy and angular momentum of the 
particle rotating on the marginally stable orbit (figure 1). 

By assumptions 1-4, the emissivity, j),, in the inertial frame of the blob can be expressed as 





ju{R) = jo exp 


(5) 


where Ris a distance from the center of the blob, and i?biob> jo> and vq are numerical constants. 
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Fig.l. Schematic picture explaining the motion ot the gas blob and the observer's plane. The xobs-axis is parallel to the equatorial plane ot the black hole 
(®), and the i/obs-axis is perpendicular to the aiobs-axis. 
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Fig. 2. Schematic picture of explaining the relationship between the blob model and the ring model. 
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2.3 Gas ring model 


Next we postulate the situation that the gas blob ereated in the innermost part of the disk is tidally 

disrupted and gets elongated by the rotational veloeity differenees in the azimuthal direetion to form 

a rotating ring (figure 2). We assume that the gas ring has the following properties: 

1. Its eross seetion has a eireular shape and its eharaeteristie thiekness is Rnng, whieh has no time 
dependenee. 

2. We negleet self-gravity of the ring. 

3. The gas ring emits radiation with monochromatie frequency, z/q, in its inertial frame. Its emissivity 
decreases with an increase of a distance from the center of the cross section of the ring (hereafter 
referred to as the center of the ring), obeying the Gaussian function. 

4. We neglect self absorption within the ring. 

5. The initial position of the center of the ring is at r = 0.98rms- 

6. The motion of the center of the ring follows the one particle orbit given by equations (3) and (4). 
The ring falls to the black hole on the equatorial plane, keeping the constant energy, E^s, and 
angular momentum, L^s- 


By assumptions 1-4, the emissivity, ji,, in the inertial frame of the ring can be expressed as 

2 " 

(5(z2-z/o), (6) 


3u{R) = jo exp 


R 


^ -^ring, 


where i? is a distance from the center of the ring, and Rrmg, jo, and z/q are numerical constants. We 
consider three cases for the ring thickness: i?ring = O.IM, 0.2M, and 0.3M. 


2.4 Methods of numerical calculations 


In order to obtain the images of infalling a blob or a ring seen by a distant observer, we solve the 
photon trajectories by the ray-tracing method (see, e.g. Karas et al. 1992). The radiative transfer 
equation is written by 


dl = g^ 



de, 


(7) 


uo J 

where / is the intensity of a ray reaching the observer, 

ju is the emissivity profile [equations (5) and (6)], di is an infinitesimal spatial interval of a 
ray in Boyer-Lindquist coordinates, and the energy-shift factor, g, is expressed as 


1 

u^{l-nA) + Wu'^' 

with 


(8) 
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w 



aAy — A [(A — ay + q] 


(9) 


A 


where and Vl{= vI"/u^) are 4-veloeity and angular veloeity of the eenter of the blob or ring, and 
A = L/E and q = Q/E are angular momenta with respect to the 0 and 9 directions per unit energy, 
respectively. Note that A and q satisfy 


A = -XobsSini, 
q = cos^ i + \^ cot^ A 


( 10 ) 

( 11 ) 


along a ray which reaches a point (xobs,|/obs) on the observer’s plane (Cunningham & Bardeen 1973). 
Here, Xobs and yohs are Cartesian coordinates on the observer’s plane, where the Xobs-axis is parallel 
to the equatorial plane of the black hole and the |/obs-axis is perpendicular to the Xobs-axis (see figure 
1). Moreover, by using polar coordinates, we take (xobs,|/obs) = (robs cos 0obs,?^obs sin 0obs), and robs 
and 0obs satisfy 0 < robs < 12M and 0 < 0obs < 27r[rad]. 

We divide the observer’s plane into x cells. Here, the number of cells is rir x = 
200 X 400, and spacings between each cell are drobs = 0.06M and robsd(/>obs = robs x vr/200. The 
distance between the center of the observer’s plane and the black hole is set to be = lOOOM. 

More detailed numerical procedures are as follows: 

1. We calculate the ray trajectories which leave each cell of the observer’s plane in the perpendicular 

direction to reach the vicinity of the black hole by applying the symplectic method (Yoshida 1993) 
obeying equations (3) and (4). Here the symplectic method is the numerical method to calculate 
test particle trajectories by solving Hamilton’s canonical equation. We assume that a photon which 
enters within the radius r = rh + is captured by the black hole and thus cannot escape from 

there. Here, rh is the radius of the event horizon of the black hole. 

2. Next, we calculate a trajectory of a test particle, which falls to a rotating black hole from r = 
0.98rjns by the Runge-Kutta method. This is to describe a motion of a center of the gas blob or 
ring. We assume that a particle which enters within the radius r = rh + 10“^M is captured by the 
black hole and thus cannot escape from there. 

3. We choose one ray, which reaches a certain cell of the observer’s plane at an observational time, 
fobs- Any photons which reach the observer’s plane at fobs after traveling time of ftravei along this 
ray was radiated by the gas blob or ring at fem = fobs — ftravei- At tem, the position and the 4-velocity 
of the center of the gas blob or ring are given by procedure 2. The emitted position of the photon 
which reaches the cell at fobs can be uniquely determined as a function of ftravei- Then by inserting 
the distance between the emitted position and a position of the center of the gas blob or ring into 
equations (5)-(6) we obtain the emissivity, jj., for each value of ftravei- 
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4. We assume the existenee of an accretion disk, which is geometrically thin and optically thick, 
outside Tins- That is, any photons that cross a disk surface should be absorbed by the disk and do 
not reach a distant observer. 

5. By using jj, calculated in procedure 3 and equations (7)-(ll), we calculate the intensity of a ray 
which reaches one cell of the observer’s plane at fobs for each value of ftravei- Here, we obtain A 
and q by inserting the position of the cell into equations (10) and (11). 

6. We then sum up all contributions by photons which were emitted at different times of fobs — ftravel 
and obtain the total intensity of rays, J, that reaches one cell at fobs- 

7. We finally integrate I over the entire observer’s plane with the area, Sohs, and divide it by S'obs to 
calculate flux, /(fobs), at each observational time, fobs, 

/(fobs) = ^J IdSobs, ( 12 ) 

where dSohs is the area of the cell. We define fmax as the time when /(fobs) reaches its maximum, 
and then, the normalized flux, F(fobs), is 

F{Us) = ( 13 ) 

Furthermore, by using jj^di at a radiative position, we calculate a photon number N{t,g)5g whose 

energy-shift factor is within g g + 5g at time, fobs^ 

N{tohs,9)dg = f f g^judidSohsdg, (14) 

iIUq J J 

where h is Planck’s constant, and the integration is made to satisfy that the energy-shift factors of 
the rays reaching the observer are within a range between g g + 5g. 

Hereafter, we re-write the observational time, fobs, simply as f. 

3 Results 

3.1 Blob model: typical case 

We first show in figure 3 the results of a typical case in which we assign (a/M,i) = (0.0,85°) and 
f^biob = O.IM. The detailed light curve around the first peak shown in panel (a) exhibits two peaks 
[see arrows in panel (a)]. The former peak appears, when the blob is at the opposite side to the 
observer’s position (phase / = vr, see figure 1), due to the gravitational leasing [effect (I) in table 1]. 
The latter peak is due to the beaming occurring at / = l.Svr [effect (II) in table 1]. After the second 
peak, the flux dumps by the redshift due to the motion of the blob moving away from the observer. 
At the time of the second peak, a centroid energy-shift factor reaches its maximum [see panel (c) of 
figure 3]. These properties are the same as the those studied of a rotating blob (cf. Karas et al. 1992; 
Dovciak et al. 2004a). 
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Fig. 3. Typical results of the blob model for the case with {a/M,i) = (0.0,85°) and iJbiob = O.lM: (a) the detailed light curve around the first peak, (b) 
the entire light curve, and (c) the time variation of the centroid energy-shift, where contributions by the photons which reach the observer after crossing the 
equatorial plane are removed in this plot. 


One might say, therefore, that the light variation in the first stage of evolution looks quite 
similar to that of a rotating blob at a fixed radius. Note, however, one important differenee between 
them; the peak intensity gradually increases with time in our model [see figure 3(b)]. Such an increase 
occurs, when a blob approaches the radius of the photon circular orbit, rph, and this phenomenon is 
caused by the focusing effect around the photon circular orbit [effect (III) in table 1]. 

Let us explain its meaning of the focusing effect around rph in a simple case. By inserting 
9 = 90° and a/M = 0 in equations (3) and (4) we find the r component of the geodesic equation to be 
[see also Shapiro & Teukolsky (1983), section 12.7]: 



where A is given by equation (10). It is easy to show that the second term on the right-hand-side of 
this equation reaches its maximum at rph = 3M. That is, when a ray travels from inside rph, the value 
of dr/dX first decreases until rph where it reaches its minimum, and then increases outward. It then 
follows that the rays emitted isotropically around rph are bound around the circular orbit for a while, 
as is illustrated in figure 4, and then, reach the observer at nearly the same time. This is why the 
observed luminosity increase occurs in figure 3(b). 

Figures 3(b) and (c) clearly prove the existence of the focusing effect around rph. The flux 
reaches its maximum at fmax, and finally decreases with time, 

while the centroid energy-shift factor of the peak is nearly constant. That is, the peak flux 
increases in the first stage because of an increase of photon numbers, not because of an increase of 
mean photon energy; which provides a support to the idea that flux increase is due to the focusing of 
rays. In the second stage (after fmax), the flux decays due to the gravitational redshift [effect (V) in 
table 1] and to an increase in the fraction of photons that are captured by the black hole. 
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Photon's path reaching observational plane 



xlM 

Fig. 4. Schematic picture expiaining the focusing around the photon oircuiar orbit, i-ph (a/M = 0.9). Each line represents a ray that finally reaches the 
observer’s plane. 

3.2 Blob model: cases with various {a,i) 

Next, the results of other eases with various eombinations of spin parameters and inelination angles 
are shown in figure 5. 

We first fix i, say, i = 85°, and let us eompare the non-spin ease [a/M = 0.0; figure 5(al)], 
middle one [a/M = 0.6, figure 5(a2)] and high spin one [a/M = 0.9, figure 5(a3)]. We immediately 
notiee that the higher the spin is, the shorter beeomes the peak interval. This is because as a increases, 
Tms decreases, and, hence, the angular frequency (fl) increases, whereas the orbital period (oc 
decreases. One may think that we can easily estimate a by measuring peak intervals on the condition 
that we precisely know M. This is, however, not so feasible in practice, since it is hard to detect 
radiation from one separate, compact blob. Rather, it is more likely that assembly of blobs (or a ring) 
fall together. Therefore, we had better pay more attention to the overall shape of the light curves. If 
we see how the maximum flux of each peak changes with time in figure 5, we find that the higher 
the spin is, the more becomes the number of peaks, and the more rapidly grows the maximum flux of 
the peak. We also notice a similar tendency in the time variations of g displayed in the lower panels. 
These facts lead us to the conclusion that how the peak flux and the energy-shift factor vary with time 
can be good indicators of the black hole spin. 

Next, we fix a and change i, finding that the flux peak increase is more appreciable in high i 
systems. In fact, the first stage evolution (when the peak flux increases) is clear when i = 85°, while 
it is not when i = 20°. This is because the focusing effect around rph is more appreciable in nearly 
edge-on systems with smaller i values. Further, the higher i is, the sharper becomes each peak in the 
g variation profile (see lower panels). This can be naturally understood, since the beaming effect is 
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tjM tjM t/M 


Fig. 5. Time variations of the normaiized flux, F{t), (upper panels) and those of the centroid energy-shift factors, g(t), (lower panels) for various combinations 
of a's and i’s. The abscissa is observational time t normalized by the black hole mass M. 


more effective in higher i systems. 
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In order to extract information regarding the spin (a/M) from the overall flux variations shown 
in figure 5, we next suppose the cases in which numerous blobs fall together. We plot in figure 6 
superposition of the light curves of numerous blobs rotating with different initial orbital phases for 
the case of a/M = 0.0 (panel al), 0.6 (a2), and 0.9 (a3), respectively. Here, we fix the inclination 
angle to be i = 85°. It is evident that the higher a is, the shorter becomes the time to reach the flux 
maximum. This is because the initial position of the blob (at ~ rms) is closer to the black hole with a 
higher a value. 

In every case, two stage evolution (rise and decay) is clear. Such a trend can be more clearly 
shown in panels bl - b3 in figure 6, in which envelope curves connecting each peak in panels al - 
a3 are displayed. To be more precise, the first stage in figure 6(bl) can be further divided to two 
sub-stages; that is, the flux first slowly increases (until t/M ~ 220) and then starts to rapidly increase 
until the flux maximum at t/M ~ 270, followed by a rapid decay at the end. Such sub-stages are not 
so clear in high spin cases [see figure 6(b3)]. 

Further, we wish to draw attention to the a-dependence of the shape of the flux curve around 
its maximum; the peak is rather broad when a/M is high [see panel (b3)], while it is sharp when a/M 
is low [see panel (bl)]. That is, the region where focusing is prominent widens as a/M increases. 
This can be understood by the consequence of the frame dragging around a rotating black hole, by 
which photons are focused in the rotating direction [effect (IV) in table 1]. 

In the high a case, the frame dragging effect plays an important role for flux peaks in panels 
(b). The each radius at which point the each flux peak is r/M = (3,2.9,2.1) = (1,1.32,1.34)rph, for 
a/M = (0,0.6,0.9), respectively. In the case of a/M = 0, the emitting radius for the flux peak is equal 
to Tph. In the case of a/M = 0.6 and 0.9, on the other hand, the radii of the flux peaks are outside of 
Tph. This can be understood by the consequence of the frame dragging around a rotating black hole, 
by which photons are focused in the rotating direction. 

We should make cautions here, however, since what are observed are not the envelope curve of 
the peaks but the total radiation energy received by an observer per unit time. We thus add three more 
panels as figures 6(cl) - (c3) displaying the time evolution of “peak area”, time integral of radiation 
energy of the peak. By comparing the bottom panels (cl) - (c3) with the middle panels (bl) - (b3), we 
understand that the sharp rise found in the middle panels are totally missing in the bottom panels and 
that such a distinction is clearer in lower spin cases. Since the bottom panels show the time integration 
of flux of each peak, while the middle panels show its maximum flux only, differences between them 
should mean differences in time width (i.e., duration) of each peak. We actually find that the time 
width of the peak significantly decreases with time because of the Doppler beaming effect is more 
enhanced when gas moves nearer to the black hole. The bottom panels will be used for our method 
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aM=0.6 
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Fig. 6. Flux variations of numerous blobs for cases with different spin parameters, a/M = o.o, 0.6 and 0.9. Panels (a) show superposition of light curves of 
numerous blobs rotating with different initial orbital phases. Panels (b) plot envelope curves connecting the peaks shown in panels (a). Panels (c) depict time 
developments of the peak areas, time integral of the peaks shown in panels (a). Here we fix i = 85°. The abscissa is observational time t normalized by the 
black hole mass M, and the ordinate is the flux normalized by its maximum value in each panel. 

of spin measurements. 
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Fig. 7. Light curves of the ring for the spins of a/M = 0.0 (left), 0.6 (middle), and 0.9 (right) and the inclination angles of i = 85° (upper), 60° (middle), and 
20° (lower), respectively. The abscissa is observational time normalized by M, and the ordinate is the flux normalized by its maximum value, F{t). We fix 

Tiring = O.IM. 

3.3 Results of ring model: light curves 

We show the light eurves of the ring for various eombinations of a and i for a fixed value of i?ring = 
O.IM in figure 7. The eontinuum eurves are nothing but the superpositions of the light eurves of 
numerous blobs with different initial orbital phases [cf. panels (e) of figure 6]. 

There are several noteworthy features found in these light eurves. First, the light eurves have 
no peaks unlike the ease of the blob model. This is beeause a distant observer at any time reeeives 
the ray whieh is at the orbital phase of 0 ~ vr (at whieh the gravitational leasing is most effeetive) or 
the one at 0 ~ l.Svr (at whieh the beaming is most effective). Next, the flux first increases due to the 
focusing effect around rph (in the first stage) and then decays due to the gravitational redshift and to 
the capture of rays by the black hole (in the second stage). These features are the same as those of the 
numerous blobs shown in figure 6. 

Let us examine the (a, i)-dependence of the light curves. Let us first fix the inclination angle 
(i) and change the spin from a/M = 0.0toa/M = 0.9 to see what changes do occur. 


15 































Fig. 8. Time-dependent properties of the radiation emitted from a ring. We fix (a/M,i) = (0.0,85°) and flri„g = o.lM. Panei (a) shows the time variations 
of the photon numbers as a function of the energy-shift factor, g. The color represents a logarithm of the photon numbers normaiized by their maximum value. 
Panel (b) shows the photon number distribution against g at the fixed time of t/M = 50. Two peaks are prominent, which are indicated by the two arrows. 
We cali ffbiue (or gred) as the g-vaiue of the blue-shifted (red-shifted) component. Panel (c) shows the time variations of guiue and gred. 


1. The higher a is, the shorter is the flux variation for all the values of i. This is beeause a ring falls 
to the blaek hole from smaller radii (reeall that the initial position at r is smaller for a higher 
spin). 

2. The higher a is, the steeper is the flux eurve in the first stage. This feature is prominent for high i 
oases; e.g., i = 85° and 60°, Note that flux inorease itself is not so appreoiable when i is low. This 
is beeause the foeusing effeot around rph, whioh is responsible for the flux inorease, is stronger 
when a is high and i is higher. 

Next, we fix a/M and ehange the inelination angle from i = 20° to i = 85°. As stated above, 
the flux inoreases more rapidly, in higher i oases beeause of more enhanoed the foousing effeot around 
Tph- Same feature is seen in other spin oases. 

3.4 Results of ring model: photon numbers 

Although the flux variation oontains useful information regarding the spin and the inelination angle, 
it is not suffioient to measure the both quantities. Another distinot information is needed. Sinee 
two pieoes of independent information are mixed in radiation flux, F{t): photon energy and photon 
numbers, it may be useful to separate them. We next investigate the (a,i)-dependenoe of the energy- 
shift faetor, g. 

Time variation of the photon number is displayed in figure 8(a). Here, the eolor represents a 
logarithm of the photon number normalized by its maximum value. We here fix {a/M,i) = (0.0,85°) 
and Rying = O.lM. Panel (b) shows the photon number profile (against g) at a fixed time of t/M = 50. 
Two peaks are prominent: one is blue-shifted peak at g ~ 1.4 and another is red-shifted peak at 
g ~ 0.5 (see arrows). Hereafter, we eall ^ibiue (or ^fred) as the g-value of the blue-shifted (red-shifted) 

16 











component. 


To consider the reason for the two peaks, let us see again the time variation of the centroid 
energy-shift of an infalling blob shown in figure 3(c). This panels exhibit two key features (or phases): 
(i) sharp peaks with large ^f-values (greater than 1.0), which occur due to the beaming, and (ii) long 
‘quiescence’ with low ^f-values (typically g ~ 0.5), which occurs when the blob is moving away from 
the observer due to the Doppler effect and to the gravitational redshift. Now we understand that 
5'biue (or g^ed) in figure 8(b) corresponds to the centroid energy-shift factor of the peaks (or of the 
quiescence) in figure 3(c). Here, we wish to stress that figure 8(b) shows a local maximum at ~ g^^d^ 
although the photon energy is low. This means that the area under the F(t) curve of the red-shifted 
component is larger than those of other energy-shift factors. 

The double-horn profile displayed in figure 8(b) resembles that of the ‘disk line’ Fe Ka emis¬ 
sion line formed on an relativistic accretion disk around a black hole (Kojima 1991; Laor 1991). We 
wish to point a big distinction, however, that our red peak is very narrow, while the red component of 
the disk line is usually broad, especially when a/M is large. This is because we are concerned with 
emission from a narrow ring (faring ^ r), while the disk line is produced on a disk with a large area. 
Superposition of emission lines from different radii produces the broad red wing of the disk line. 

Finally the time variations of these quantities are shown in Panel (c). The spin and inclination- 
angle dependence is expected, which will be shown in next subsection. 


3.5 Results of ring model: energy-shift factor 

We show the time variations of ^^biue and ^fred for various combinations of a and i in figure 9. First, 
we immediately notice a clear trend that ^f-values depend more sensitively on i, rather than a, and that 
the higher i is, the longer becomes the energy-shift interval between the two lines. This can be easily 
understood, since the Doppler is more pronounced when i is higher. Let us next fix i and see how the 
results change by decreasing a from high a/M{= 0.9) values to lower one. Here is a summary: (i) 
When i is low; e.g., i = 20°, both of r^biue and ^^red decrease with a at low i, since gravitational redshift 
is stronger in order that the initial position of the ring (f» r^s) is smaller [see figures 9(cl) and (c3)]. 
(ii) When i is large; say, e.g., i = 85°, guue does not depend so much on a, while since the beaming 
effect is stronger than the gravitational redshift, while ^fred is smaller due to the gravitational redshift 
and the redshift due to the motion of the ring fraction moving away from the observer. 
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Fig. 9. Time variations of gbiue (dot line) and gred (solid line) for a variety of o/M = o.o (left), 0.6 (middle), and 0.9 (right) and i = 85° (upper) ,60° (middle), 
and 20° (bottom), respectively. Each curve is plotted until the time when the normalized flux decays to reach F = 1/2 (i.e., half of the maximum flux). 


4 New method of black-hole spin measurement 

4.1 Introduction of new quantities 


As is demonstrated in the previous seetion, radiation from an infalling gas ring contains useful infor¬ 
mation for the a and i measurements of a central black hole. In this section we propose a methodology 
for the measurement of a and i. In the proposed method we use the following three quantities: (i) 
gredit), energy-shift factor of the red-shifted spectral component, (ii) 5'biue(^), the same but of the 
blue-shifted component, and (iii) F{t), the normalized radiation flux (see figures 7 and 9). From 
these quantities we calculate the following six quantities: 


2 /*^inax 

r 

3M Jt„ 

A = Al -f A2, 


m - m) 

1 - F{ti) 


dt, 


F{t)dt, 


(15) 

(16) 
(17) 
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(18) 


B 

C 

D 


II- 9redit)dt 

{t{ - ti) 
fti dhluei.^'jdt 
{tf - ti) 


(S'red)) 


(S'blue)) 


(19) 

( 20 ) 


where we fix the initial time to be ti/M = 60 and tf is defined in sueh a way that F(tf) = 0.5. Note 
also that t is normalized by M. One may think that this definition of ti looks quite arbitrary and that 
the value of Ai erucially depends on how to ehoose ti. Fortunately, however, this is not the ease, sinee 
only the flux which is substantially greater than F{ti) contributes to the integral in Ai. In other words, 
Ai depends not critically on ti but on t[, when the flux starts to grow rapidly. 

Let us next consider the physical meaning of Ai, A 2 and A by taking an example of a simple 
functional form for F{t). Namely, we assume that F{t) stays nearly constant F{t) ~ F{ti) until t[ 
and that it grows as a linear function of t until fmax; 

F(t)»F(ii) + |l-F(ii)l^, (21) 

^rise 

where frise = f^max — 1[ is the typical rise time. We then find 


Ai ^ frise/M = (t„,, - f') /M. (22) 

Note that Ai does not depend on the initial flux, F{ti) nor ti. If we assume the linear functional form 
for F(t) also for a decay, 

F{t)^F{t,) + [l-Fit,)]^f^, (23) 

^decay 

where fdecay = t{ — fmax IS the typical decay time, we obtain 


A2 ^ (ff-fmax)/iVf 


We finally have 


(24) 


2l = 2li + 2l2 ^ {ti-t’^/M 

Thus the value of A means the typical timescale of flux variation. 

In table 2 we summarize the physical meanings of the quantities, A-D. 


(25) 


In the next section, we calculate Ai-D for a variety of combinations of a, i, and Rnng to seek 
for the best combinations of the quantities to extract information on the spin from observational data. 
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Table 2. Physical meanings of the quantities, A-D. 


quantity 

physical meaning 

equation number 

A 

timescale of flux variation 

(17) and (25) 

B 

time averaged gfi-ed it) 

(18) 

C 

time averaged 5rblue(f) 

(19) 

D 

relative difference between B and C 

(20) 



0.3 0.6 0.9 1.2 1.5 

C 


Fig. 10. Relation between A and C for the inclination angles of i=20, 40, 60, 80 and 85° and for the spin parameters of a/M = 
-0.9, -0.8, -0.6, -0.4, -0.2,0,0.2,0.4,0.6,0.8 and 0.9. We fixed the ring width to be Firing = O.lM. The data points with the same a value are 
connected by the lines. The two arrows indicate the direction of increasing a. 




Ai B C 

Fig. 11. Relation between the quantities except for A and C for the same inclination angles, spin and ring width on figure 10. 
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4.2 Dependence of Ai-D on (a,i) 

The calculation results are summarized in figure 10 and 11. It is important to note that all the curves 
in each panel do not cross each other, that is; we can, in principle, uniquely determine a and i, once 
the values of Ai-D are specified. 

Particularly, we notice that the lines (on which a is constant) are nearly horizontal; that is 
A does not so much depend on i (figure 10). This fact leads to an important conclusion that we can 
roughly estimate a solely by A. We emphasize that it shows a way to break the classic a-i degeneracies 
of the spectral methods of measuring spin, and is not particularly model dependent like most QPO 
methods. 

Let a/M increase from 0.0 to 0.9 for a fixed value of i to see how the key quantities of Ai-D 

vary. 

1. We find that all of the A values decrease with an increase of a/M. This is because the timescale 
of flux variation decreases with an increase of a/M both in the rise and decay phase for all i (see 
figure 7). 

2. With an increase of a/M both of B{= (gred)) and C{= (ghiue)) decrease, unless i is very high, 
because of the gravitational redshift (see figure 9). Since the decrease of C is slower than that 
of B, D[= {C/B) — 1] increases with an increase of a/M. When i is very high (e.g., i = 85°), 
conversely, C rather increases (while B still decreases) with an increase of a/M, This is due to the 
beaming effect, which is more effective for edge-on cases 

4.3 Dependence of Ai-D on the thickness of the ring 

In the previous subsection, we propose a new method for the spin measurement. One may think, 
however, that the results may sensitively depend on the ring shape, such as the characteristic thickness 
of the ring. To demonstrate that this is not the case, we need to consider how the values of Ai-D 
depend on the thickness of the ring. We will demonstrate that the results will not depend on Turing as 
long as if we use A, instead of Ai or A 2 . 

We display in figure 12 the relation between several quantities for i?ring = 0.2M [panels (al) 
and (a2)] and 0.3M [panels (bl) and (b2)], respectively. By comparing these panels with figures 11(a) 
and 10, we understand that the relation between A and C do not critically depend on the thickness of 
the ring. Since we put the center of the ring (where the emissivity ji, is maximum) at the same position 
for all the models, we can easily understand that B-D do not so much depend on a thickness of the 
ring. For low inclination angles (say, i = 20°), hence, we can determine a and i by using B,C and D. 
For high inclination angles (say, i > 60°), however, the photon number of the red-shifted component 
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is much less than that of the blue-shifted one so that the precise measurements of from aetual 
observations might be diffieult. From heneeforth we mainly use the information A and C. 

Figures 12(al) and (bl) display how the relation between Ai and A 2 depend on the thiekness 
of the ring ehanges. Here, we fixed a and i. Comparing these with figure 9(a) showing the ease with 
-Rring = O.IM, we see a clear tendency that the larger -Rring is, the larger is Ai and the smaller is 
A 2 . This ean be understood, sinee a larger thiekness means a longer emission timeseale, thereby Ai 
being increased. (Note that the motion of the center of the ring is kept the same regardless of the ring 
thiekness, but the epoeh when the last tip of the ring material passes a eertain radius (say, rph) is more 
and more delayed, as the ring thiekness inereases.) Since the deeay timescale deereases, so does A 2 , 
as i?ring inereases. 

To summarize, the thicker the ring is, the larger Ai is and the smaller A 2 is, whereas A (= 
Ai + A 2 ) does not eritieally depend on the ring thiekness [see figure 10, figures 12(bl) and (b2)]. 
In fact, the relation between A and C shown in figures 12(al), (bl) and figure 11(a) exhibit similar 
tendeneies among different ring models. 

Let us finally demonstrate that we can certainly have good guesses for a/M and i, by using 
the {A, C) diagram in figure 12. The results are summarized in table 3. We find that the differenees 
in a/M and i among different models are within ranges of ~0.01 — 0.1 and ~ 1°, respectively. To 
eonelude, we ean measure the spin and inelination angle with good aeouraey by using A and C. 
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(i) iJring = 0.2J/ 




(ii) -Rriiig = 0.3,1/ 



Fig. 12. Panels (a1) and (b1) [or panels (a2) and (b2)] are the same as figure 11(a) [or figure 10] but for different values the ring thickness: Rring = 0.2M 
[(a1) and (a2)] and 0.3M [(b1) and (b2)], respectively. 
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Table 3. Estimated values of a and i from given A and C. 


quantities 

thickness of the ring 

spin 

inclination angle 

{A,C) 

Rring/M 

afM 

i° 

(100,0.9) 

0.1 

0.72 

45 


0.2 

0.74 

46 


0.3 

0.76 

46 

(100,1.2) 

0.1 

0.78 

69 


0.2 

0.77 

70 


0.3 

0.76 

69 

(300,0.9) 

0.1 

-0.44 

22 


0.2 

-0.43 

22 


0.3 

-0.44 

23 

(300,1.2) 

0.1 

-0.27 

64 


0.2 

-0.27 

64 


0.3 

-0.28 

65 


5 Discussion 

5.1 Brief summary 

In this paper, we propose a new method for measuring a and i from observational data for given 
mass M. In this method, we consider a gas cloud, which is assumed to have a ring shape, falling 
onto a black hole. Assuming that the gas ring emits monochromatic radiation, we calculated the flux 
variation and the time averaged frequency shift of the radiation received by a distant observer for 
various combinations of a and i. We have demonstrated that we can uniquely determine a and i from 
observational data, not critically depending on the thickness of the ring. 

The actual procedures of our method are in three steps: 

1. Obtain time-sequence data of normalized flux F{t) and photon numbers N{t,g)5g as functions of 
frequency or g (frequency-shift factor). 

2. Calculate the key quantities, A and C, by equations (15)-(17) and (19). 

3. Estimate a and i by using figure 10. 

In addition , we can also estimate the ring thickness from Ai and A 2 by using figure 11(a), figures 
12(al) and (bl). Even when the frequency information is not available, we still have a good guess on 
a solely from A, especially when the ring thickness is relatively large (~ 0.3M). 
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5.2 Distinctive features of the proposed method 

The most distinctive feature of the proposed method is that we can directly measure a black hole 
spin. This is made possible by focusing on non-periodic light variation events in the innermost region 
at r < rjns. By contrast, the widely used methods are based on the steady or periodic phenomena 
occurring at r > rms- In this respect our method is not only independent and complementary to other 
methods, but also advantageous in probing the innermost space-time structure around a black hole. 

The most popular methods of spin measurements based on continuum or Fe Ka line spectra 
totally depend on the assumption that the inner edge of the disk is at r^s (Duro et al. 2011; McClintock 
et al. 2011; Steiner et al. 2011; McClintock et al. 2014; Reynolds 2014). We wish to point that there 
is no solid physical reason for this assumption; i.e., the disk could be truncated at a radius greater 
than Tins- It is of great importance to note the facts that line and continuum spectra are sensitive more 
to the radius of the inner edge of the emission region than to the spin values. This is obvious for the 
continuum spectra, since the smaller rms is, the higher is the maximum disk temperature (Tmax) and 
the smaller is the region emitting blackbody radiation with while the shape of the gravitational 
potential at r > r^s does not largely change with a change in a. This is also true for line spectra. 
Kojima (1991) compared two line profiles, one from a disk with the inner edge at 6M around a 
non-rotating black hole and another from a disk with the same radius of the inner edge but around a 
maximally rotating black hole, finding little difference between them (see his figure 2). In short, these 
major methods only indirectly measure spins. 

One may think that our method might crucially depend how to take the initial time (ti) and 
radius, when and where a rotating ring starts to fall. This is not the case, however. We can arbitrarily 
choose the initial radius as long as r as we have argued in subsection 4.1. This is because the 

gas ring is initially rotating on nearly a circular orbit and so its emitted flux is nearly constant at first, 
not contributing to Ai. As a result, A does not critically depend on the initial position nor time. 

5.3 Observational implications: shot analysis 

Our ring model may possibly be applied to the X-ray shots of black hole binaries. X-ray shots are 
flare-like light variations with sharp peaks (see Oda et al. 1971; Gierlifiski & Zdziarski 2003; Warren 
et al. 2005) and are observed during the so called low/hard state, which is characterized by power-law 
spectra and is associated with substantial variations (see e.g., a review by Done et al. 2007). 

Negoro et al. (1994) developed the technique of superposed shots, adding plenty of shot 
profiles (in X-ray bands) by aligning their peaks (called ’’superposed shot” analysis: see also Liu & Li 
2004; Wu et al. 2007; Yamada et al. 2013). Furthermore Feng et al. (1999) detected the shot profile 
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with high time resolution, and found that (1) the time profile of the superposed shot is rather time 
symmetrie with respeet to its peak time and that (2) it ean well be fit with the sum of two exponential 
functions: one with time constants of ~ 0.01 [s] and another of ~ 0.1 [s]. 

The expected variation timescale by the ring model is 

tvar ~ 1O(M/1OM0) ms (26) 

or ~ 0.01 s for Cyg X-1 with 12.8 < M/Mq < 14.8 (Orosz et al 2011). This means, the light 
variations predicted by our ring model may be contained within the superposed shot. The timescale 
is more favorably long for supermassive black holes. Therefore, it might be possible to determine a 
and i by future observations. 

By means of global MHD simulations Machida et al. (2003) studied the time evolution of a 
torus, which is initially threaded by a weak toroidal magnetic field, around a non-rotating black hole. 
They carefully analyzed the simulation results and found that gas clouds (or streams) with a spiral 
shape intermittently fall onto the black hole from the inner edge of the torus. These gas clouds are 
heated by dissipation of magnetic energy via magnetic reconnection and are thus expected to emit 
intense X-rays. This is exactly the situation that we have postulated in our blob or ring model. 

We, however, note that the shape of accreting gas clouds is not a full ring but is a spiral 
according to the MHD simulation. It is thus necessary to check whether our proposed model based 
on the full rings can apply to spiral-shaped gas clouds or not. We calculate a half ring model as the 
intermediate case between the blob model and the ring model and show its light curves in figure 13. 
The light curves resemble those of the blob model in the sense that they have several peaks (or humps) 
due to the gravitational leasing and the beaming, but the duration of each peak is much longer in the 
half-ring model. We also confirm the tendency that the peak flux of the half ring first increases by the 
focusing effect around rph and then decays due to gravitational redshift. Such overall variations are 
similar to those in the ring model. We thus expect that the superposed light curves of the half rings 
(with different initial phases) are equivalent to those of a full ring (recall figure 6). In conclusion, we 
can use the ring model for spin measurements regardless of the precise shape of the gas cloud; it can 
be a blob, a half ring, or a full ring. 

5.4 Observational implications: feasibility 

How accurately can we measure the spin value when applying our method to the actual X-ray ob¬ 
servational data? In subsection 4.3, we have demonstrated that we can estimate the spin value by 
measuring A. In this subsection, we evaluate an accuracy of the spin measurement by performing 
Monte Carlo simulation for light curve of the gas ring. 
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a/AI = 0 a/M = 0.9 



Fig. 13. Light curves of a half ring (by the solid line) and those of a full ring (by the dashed line) for a/M = 0 (left panels) and 0.9 (right panels), and for i=85° 
(upper) and 20° (lower), respectively. We assumed R^ns = O.lM in all models. 


Let us first examine the case of observing Cyg X-1. The time resolution of the RXTE is 1 [ms] 
(Feng et al. 1999; Liu & Li 2004; Wu et al. 2007). This timescale corresponds to 0.07tvar for the 
black hole mass of M = 14Mq (Orosz et al 2011), where, fvar is the time scale of flux variation [see 
equation (26)]. Hereafter we thus assign a constant time bin; At = O.OTtvar- Suppose that we observe 
Cyg X-1 with the RXTE satellite. According to Feng et al. (1999), however, the number of photons 
that are received during the time bin of O.OTtvar is only 2 [counts/keV] even at the peak. We thus give 
up using the RXTE data and calculate the cases of observations with future satellites, such as Athena 
(or EOFT). Since the effective area of these missions is about 10 (100) times as large as that of the 
RXTE, the average count rate around the peak is = 20 (200) [counts/(0.07fvar)/keV]. 

The actual procedures of evaluating the spin parameter ranges (Ogst) are as follows: 

1. We first choose one value of the spin parameter (otrue) and, assuming that the photon number at the 
flux maximum to be A^ax given above, we calculate the expected (averaged) number of photons 
[A(f)] received per each time bin with a length of 0.07fvar, noting N{t) (x F{t). 

2. On the averaged photon number profile [A(t)] we superpose random errors [5Aj(f)] (with i = 
1,2, • • •), assuming that photon statistics obeys Poisson distribution [with the standard deviation of 
o'N = \/N(t)]. Figure 14 give one such example. 

3. We superpose 500 such shot profiles. 
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Fig. 14. An example of the time profile of the observational flux, in which random errors are superposed by assuming the Poisson statistics. The solid curve 
represents F(t), the result of figure 7, while the vertical bars represent one example of N{t) + SNi{t). Here, we postulate the case that we observe Cyg 
X-1 with Athena, and the model parameters are (otrue/M, i) = (0,85°). The time interval corresponds to 0.07tvar — 14M. 


Table 4. estimated spin values of Cyg X-1. 


Cttrue/ M 

Athena 

flest fori = 85[°]) 

i = 40[°]) 

LOFT 

* = 85[°]) 

i = 40[°]) 

0.8 

0.77-0.83 

0.68-0.85 

0.78-0.81 

0.79-0.85 

0.4 

0.31-0.46 

0.28-0.50 

0.35-0.40 

0.38-0.50 

0 

-0.15-0.11 

-0.15-0.10 

-0.10-0.02 

-0.01-0.09 

-0.4 

-0.54- -0.22 

-0.51--0.31 

-0.48- -0.35 

-0.49- -0.35 

-0.8 

-1.0--0.61 

-1.0 -0.70 

-0.87- -0.74 

-0.88- -0.74 


4. We then normalize the time sequence of the number counts [Ntot (^)] of each run by their maximum 
value, and call it the re-defined flux distribution, F{t). 

5. We calculate A by inserting F{t) to F{t) in equations (15)-(17) and estimate the spin parameter 
(Oest) by using figure 10. 

6. We repeat the same process for 10^ times, plot the distribution of the estimated spin, Ogst, and 
evaluate its l-a confidence range. 

We summarize the results in table 4 for each of observations with Athena and LOFT, re¬ 
spectively. In the case of Athena (or LOFT), we can determine the spin value with following error, 
5a/M = 0.06 - 0.39 (0.03 - 0.14). 

Similar analysis of the case of Cyg X-1, we estimate the accuracy of spin in the case of MCG- 
6-30-15. The time resolution of RXTE is 10^[s]= O.lfvar, where we postulate the mass is lO^M© (Lee 
et al. 2000; Reynolds 2000). Then we set to At = O.ltvar- The averaged photon count rate of one shot 
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Table 5. estimated spin values of MCG-6-30-15. 



RXTE 

aest for i = 85[°]) 

i = 40[°]) 

Athena 

* = 85[°]) 

i = 40[°]) 

LOFT 

i = 85[°]) 

i = 40[°]) 

0.8 

0.74-0.87 

0.77-0.83 

0.68-0.85 

0.73-0.80 

0.78-0.81 

0.79-0.85 

0.4 

0.21-0.58 

0.32-0.46 

0.28-0.50 

0.34-0.49 

0.35-0.40 

0.38-0.50 

0 

-0.30-0.33 

-0.13-0.10 

-0.16-0.10 

-0.15-0.13 

-0.1-0.02 

-0.01-0.09 

-0.4 

-0.70-0.08 

-0.53- -0.16 

-0.55- -0.22 

-0.51- -0.31 

-0.48- -0.34 

-0.49- -0.35 

-0.8 

-1.0--0.25 

-1 - -0.56 

-1.0--0.55 

-1.0--0.70 

-0.88- -0.74 

-0.88- -0.74 


peak is 100 [counts/(0.1tvar)/keV]. Further, we set the number of superposed shots to 10[shot], sinee 
the time interval between start and end time of the observation, 10^ [s] is about 10 times as large as 
an time duration of the shots (Lee et al. 2000; Reynolds 2000). In addition, sinee the effeetive area of 
Athena (or LOFT) is about 10 (100) times as large as that of the RXTE, the average eount rate around 
the peak is Wmax = 10^ (10"^) [eounts/(0.1fvar)/keV]. 

The results for MCG-6-30-15 are summarized in table 5. By using Athena (or LOFT), we can 
estimate spin value with following error, 6a/M = 0.07 — 0.45 (5a/M = 0.03 — 0.14), while in the case 
of RXTE, we can scarcely detect the spin value with the rough accuracy, 5a/M = 0.13 — 0.78. 

5.5 Remaining issues 

We have so far argued that the proposed method should in principle work, however we recognize 
number of essential issues to be considered before applying this method to actual observational data. 
We first evaluate what resolutions are necessary for spin determination. 

The first issue is time resolution. As we indicated by equation (26), the time resolution re¬ 
quired by the ring model is less than 1O(M/1OM0) ms. This value itself is attainable but we need 
good photon statistics, which is another concern. Note that this timescale is much longer than the 
dynamical timescale (or the period of orbital motion), which is on the order of less than 1 ms. This is 
because we are considering gas falling with keeping a finite angular momentum. 

The next one is energy resolution. In order to estimate i from the G-value within an error 
of 20°, the required energy resolution is 0.1 x Enne from figure 10, where Enne is the energy of the 
fluorescent Fe Ka line in the inertial frames of the emitting gas. 

Finally, we need to know M with high accuracy to estimate a accurately. As we mentioned in 
subsection 4.2, we can measure a mainly through A, and A depends on M [see equations (15)-(17)]. 
Let us denote 5M as an error in M estimation. Then the error in A is evaluated as 
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Table 6. relationship between the estimated spin 


and true one for a given error in mass. 


mass error 

true spin 

estimated spin 

\6M\/M 


<^est /^ 

0.1 

0.8 

0.74-0.85 


0.4 

0.30 - 0.49 


0 

-0.17-0.12 


-0.4 

-0.55 - -0.23 


-0.8 

-1.0--0.62 

0.2 

0.8 

0.69 - 0.92 


0.4 

0.21-0.60 


0 

-0.31-0.28 


-0.4 

-0.68 - -0.021 


-0.8 

-1.0--0.40 


5A _5M f 5M\ ^ 


(27) 


Using this relation we ealeulate the range of the estimated spin, Oest, for a given (true) spin value, Otrue, 
for two eases with 5M/M = 0.1 and 0.2. The results are listed in table 6. We see in this table that 
we can determine a/M with good accuracy of within ~ 0.05 for high spin cases, or less than 0.20 for 
low spin cases (say, a ~ 0), if the mass determination accuracy is 5M/M = 0.1. When 5M/M = 0.2 
the situation gets worse, but the spin determination accuracy is still ~ 0.3 even for low spin cases. 

Finally, we give a summary of future issues. 


1. We considered particular shape of infalling gas cloud and assumed no changes of the shape with 
time. That is, a distortion of the ring shape in a deep potential well was not considered. We need 
to relax this assumption and improve the method by using, e.g., MHD simulation data. 

2. The temperature and the emissivity (jo) of the ring is kept constant. This may not be so realistic, 
since heating by magnetic energy release and/or radiative cooling should be effective. 

3. We assumed monochromatic radiation form the gas clouds for the sake of simplicity. This is a 
critical assumption, on which the proposed method is constructed, but its physical reason is not yet 
clear. We need to expand our methodology so as to incorporate more general cases of continuum 
radiation etc. 

4. Actual correspondence with real observational data is under consideration. It will be important to 
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identify which data really represent an infall of a gas ring in connection with real observational 
data and MHD simulation data. 
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